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^^ We show that the "ridge" phenomenon in the two-particle angular correlation function, 

y^N as observed by the CMS experiment, can be reproduced by implementing an impact pa- 

j rameter dependent azimuthal correlation of the scattering planes of individual partonic 

interactions. Such an approach is motivated by the observation that even for moderate 

i""i impact parameters a substantial number of partonic interactions may be produced, while 

^^ at the same time the protons are sufficiently far apart to create a preferential direction in 

^T^ azimuth. 

Pl^ A re-tune of the PythiaG Z2 tune based on underlying event and minimum bias distribu- 

QJ tions measured at the LHC shows that a better description of data can be obtained with 

j^ , this approach and that some tension existing between underlying event and minimum bias 

distributions can be removed. We show that, even though the CMS result on the angular 

^S| correlation function itself is not used in the re-tune, we can predict the appearance of 

^ long-range, near-side angular correlations in proton-proton collisions. 

00 

O 1 Introduction 

(T^ The CMS ridge effect is a two-particle angular correlation effect observed by the Compact Muon 

^^ Solenoid (CMS) experiment in Large Hadron Collider (LHC) proton-proton collisions at high 

^^ charged track multiplicities {Nch > HO) and in a specific transverse momentum range {pt = 

1 — 3 GeV). This effect, together with other visible structures in the R{Ar]^A(p) correlation 

^ distribution, was described in 1 1 . 



Two-particle correlation function R ( Ar^, A^) is defined as: 

Data is binned according to charged track multiplicity Nch- The signal 6'Ar(A77, A(p) consists of 
the charged two-particle density, while the background B^iAr], A(j)) is given by the distribution 
of uncorrelated particle pairs - the product of two single-particle distributions. Finally, the data 
is averaged, weighted with bin multiplicity, over all bins. The analysis is repeated for four sets of 
data. On one hand two minimum bias sets (all Nch)^ one including all particles with transverse 
momenta above 0.1 GeV and the other including all particles with transverse momenta between 
1 and 3 GeV. On the other hand two high- multiplicity sets {Nch > HO), again with the same 
two transerverse momentum selections. 
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Some of the effects reported in [P, including the near-side peak at (Ar^, A^) = (0,0), the 
away-side ridge at (A77, A^) = (At^, tt) and the Gaussian ridge at (A77, A^) = (0, A^) can be 
explained with single two-to-two partonic interactions. The first two are visible in all sets of 
data, while the third one is most clear in the minimum bias pt > 0.1 GeV case. A fourth effect, 
the near-side ridge, a long-range azimuthal correlation at (A?7, Acj)) = (A77, 0) only visible in 
high- multiplicity events at moderate pt, requires further study. It is for this last effect that we 
propose a model. 

For our study we observe the effect of our modification of the PythiaG [2^ Monte Carlo (MC) 
event generator on select observables and consider changes in a few existing PythiaG param- 
eters to counteract the side-effects of our modification. This latter step can be considered a 
re-tuning to CMS data. Note that we only use a limited set of CMS data and start from the 
existing PythiaG tune Z2 [3 . More global tuning including other experiments' data was not 
within the scope of this study, but may be added later. 

2 The azimuthal ahgnment model 

For large enough impact parameter b (figure fl]), the multiparton interactions in proton-proton 
collisions tend to lie in the collision plane of the hardest interaction and the final state particles 
will have similar azimuthal angle (/) - this results in near-side effects. Furthermore, an explana- 
tion for the ridge effect with multiparton interactions would require enough such interactions 
to be taking place, which leads to high-multiplicity events. At the same time we require that 
the multiparton interactions are semi-hard, and thus yield moderate-p^ particles. Finally, we 
are dealing with incoming partons with very different x^j and as such will have interactions in 
a broad pseudo-rapidity range r] - this gives rise to long-range effects. So far, everything is still 
consistent with the observations made by CMS. 




(transverse plane) 

Figure 1: Protons separated by impact parameter b. 

What is still a problem, is that high-multiplicity events are generally central collisions which 
have an impact parameter 6^0, while the definition of the collision plane of the hardest 
interaction requires large b. In light of this issue, we study whether a small upward fluctuation 
in the amount of multiparton interactions, for the case of moderate impact parameter, suffices 
to explain the CMS ridge effect. 

The modification we introduce goes on top of the most recent multiparton interaction model 
currently in PythiaG [4 . In this existing model, the amount of multiparton interactions, a 
measure for the activity, is inversely proportional to impact parameter b (VINT (139), rescaled 
to bavg — 1 for the minimum bias case). The azimuthal angle (\) (VINT(24)) is chosen randomly. 
This last point makes that angular correlations - also the long-range, near-side ones - would 
be missing in events generated with PythiaG. 
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We propose sampling random points (x^, yi) in Gaussian proton profiles (figure |2|, these protons 
being separated by impact parameter 6, and using trigonometry to calculate the (/)-offset from 
the hardest interaction. To allow for some tuning freedom we add a scaling parameter a to the 
impact parameter b. Ideally, the scaling parameter would be one. This results in: 

^i = (t^hardest + arctau I — r I (2) 

V(X2 +a- b/bavg) - XiJ 

We implement the modification for two different modes of the multiparton interaction model 
of PythiaG which both make use of hadronic overlap according to Gaussian distributions. In 
those cases, the above (/)- definition makes sense. In our tuning activity reported in section [3J 
we focus on the mode which uses double-gaussian matter profiles (MSTP(82) = 4). 




Figure 2: Sampling of random points in Gaussian proton profiles, separated by impact param- 
eter 6, and introduction of scaling parameter a to allow some tuning freedom. 

The modification has several implications. We study two sets of data: CMS underlying event 
(UE) data [3 , showing the charged multiplicity Nch and transverse momentum sum ^pr in 
the region transverse to a jet or hard interaction (figure [3|, and CMS minimum bias (MB) 
data [5 , showing the charged multiplicity Nch integrated over azimuthal angle (/). By introduc- 
ing the modification, we generate interactions with an azimuthal separation from the hardest 
interaction smaller than would be the case with the previous uniform azimuthal distribution. 
The interactions get shifted to the toward/away regions and the plateau for TV^j^^^^^erse (^^.Qpg 
(figure [4J top). Re-raising this plateau to describe the data requires a re-tune, modifying the 
PT-cutoff and by proxy the activity, Nch- The pT-cutoff in PythiaG is given by: 

pjp^n^EcM)=p'T-(^) =PARP{82)-(^) (3) 

\rjREF J \^REF J 

In contrast with the clear effect on UE results, we expect little or no sensitivity to the modi- 
fication for the MB results, which are integrated over azimuth (j) (figure [4J bottom). Possibly 
this diffence in sensitivity also allows to lift some of the tension which exists between the UE 
and MB descriptions. 

Leading Trackf Jet Direction 




Figure 3: Areas in A^ with respect to the leading track jet. 
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Nch at Vs = 0.9 TeV, |)/| < 2.0, 60 < \A(p\ < 120 




Nch at Vs = 7-0 TeV, \>j\ < 2.0, 60 < |A<f | < 120 



Leading Track-jet p^[Gfy/c] 



Nchat\/s = 0.9TeV, If/I < 2.4 





Figure 4: Overview of the sensitivity of Nch observables in CMS UE (top) and MB (bottom) 
data to changes in scahng parameter a [purple], with Z2 [black] as a reference. 

3 Tuning 

In this second part we report the results of our small-scale automized three-parameter re-tuning 
to the two earlier described CMS data sets. We start off with a review of the sensitivity of 
the observables to the three parameters and end with two tunes, one simple tune to just four 
Nch distributions (transverse & total Nch and 0.9 & 7.0 TeV) and one two-step tune to all MB 
and UE observables in the data sets. The first tune allows us to get a feeling of the parameter 
space in play, while the second one aims to reach a more solid result fixing the pT-cutoff based 
on the MB observables (insensitive to scaling parameter a) and then using the UE observables 
to fix a. For the actual tuning, we make use of the PROFESSOR package P, which takes care 
of automated tuning based on Rivet plots [7] of observables with reference data. For the 
interpretation of the tune result we make the comparison with existing tune Z2* rather than 
with tune Z2. Z2* is a Professor re-tuning of Z2 to CMS UE data, for parameters PARP(82) 
and PARP(90). 

3.1 Sensitivity 

The first observation (figure pi) is that PythiaG p^-reference PARP(82) and energy-scaling 
parameter PARP(90) affect all activity, both transverse and total activity at both 0.9 and 
7.0 TeV. For variations in the p^-reference (red/blue solid) the effect is the same at both 
energies, while for variations in the energy-scaling parameter (green/orange dashed) the effect 
is opposite at the two energies. This is to be expected since 0.9 and 7.0 TeV lie on both sides 
of reference energy 1.8 TeV used in the pT-cutoff formula. On the part of a (figure ffl, we can 
see a clear effect in the transverse region (UE dataset) and little to no effect in the ^-integrated 
case (MB dataset). 
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Nch at Vs = 0.9 TeV, \i]\ < 2.0, 60 < \A(p\ < 120, p!p^^ > 3GeV/c 

' ' ' I ' 

CMS UE data 

Z2 (1.83, 0.28, 

Z2 + low PARP(82) (1.70, 0.28, 0.00) 

Z2 + high PARP(82) (1.95, 0.28, 0.00) 

Z2 + low PARP(9o) (1.83, 0.20, 0.00) 

Z2 + high PARP(9o) (1.83, 0.35, 0.00) 

, , , ^ 



Nchat v/i=0.9TeV, |^| < 2.4 



Nch at ^/s = 7.0 TeV, \tj\ < 2.0, 60 < \A(p\ < 120, pf""" > 3GeV/c 





F 




- CMS MB data ^' ^---j 

- Z2 (1.83, 0.28, 0.00) b^^ , 

- Z2 + low PARP(82) (1.70, 0.28, 0.00) L_|__ 

- Z2 + high PARP(82) (1.95, 0.28, 0.00) L. 

- Z2 + low PARP(9o) (1.83, 0.20, 0.00) 

- Z2 + high PARP(9o) (1.83, 0.35, 0.00) 





- CMS UE data 

- Z2 (1.83, 0.28, 0.00) 

- Z2 + low PARP(82) (1.70, 0.28, 0.00) 

- Z2 + high PARP(82) (1.95, 0.28, 0.00) 

- Z2 + low PARP(9o) (1.83, 0.20, 0.00) 

- Z2 + high PARP(9o) (1.83, 0.35, 0.00) 



Lf^ 




Nch at \/J = 7.0 TeV, |i?| < 2.4 




- CMS MB data 

- Z2 (1.83, 0.28, 0.00) 

- Z2 + low PARP(82) (1.70, 0.28, 0.00) 

- Z2 + high PARP(82) (1.95, 0.28, 0.00) 

- Z2 + low PARP(9o) (1.83, 0.20, 0.00) 

- Z2 + high PARP(9o) (1.83, 0.35, 0.00) 




Figure 5: Overview of the sensitivity of Nch observables in CMS UE (top) and MB (bottom) 
data to changes in PARP(82) [blue, red sohd] and PARP(90) [green, orange dashed], with Z2 
[black] as a reference. In brackets the values (PARP(82), PARP(90), a) are marked. 



3.2 Tune one 

This first crude tune, which we call Z2R, is made with the PROFESSOR package to just four 
observables (two UE and two MB) and confirms what is also more or less visible by eye. We 
want a moderate a, so as not to destroy the match with data so much that it cannot be restored 
(either in MB, UE or both), while still having enough power to introduce the intended long- 
range near-side effect. Next, we need a slight lowering in the p^-reference (compared to optimal 
tune Z2*) to re-raise the Nch plateau in the transverse region. The energy-dependence will be 
of less importance. We find exactly this in our Z2R Professor tune (table fl]), for which we 
used the cubic interpolation mode. In general we find that, to begin with, the match with 
data for the four observables to which we tuned, is of the same quality as in the case of Z2*. 
Furthermore, also for the observables in the dataset which we did not include in the tuning, 
the match remains acceptable. We show the graphical result of tune Z2R (red solid) in figure 
pi with tunes Z2 (black dotted) and Z2* (blue dashed) as reference. 





PARP(82) 


PARP(90) 


a 


Z2 


1.83 


0.28 


0.00 


Z2* 


1.93 


0.23 


0.00 


Z2R 


1.87 


0.23 


4.15 



Table 1: Result of the 1-step 3-parameter tune to 4 observables. 
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3.3 Tune two 

The second tune we consider is made in two-steps, we call it Z2R'. This time using all the 
observables in the same UE and MB CMS datasets, we again perform an automized Professor 
tune with cubic interpolation. In the first step, the pT-cutoS (both PARP(82) and PARP(90)) 
is fixed to the MB data, disregarding any match with UE data. In the second step, a is tuned 
to the UE data. After the first step, the match with data is good for MB, but less so for UE. 
After the second step, also the match with UE is restored to an acceptable level, comparable 
to the Z2* tune. Quantitatively, we again find the tune to be insensitive to PARP(90), while 
PARP(82) and a settle on values inbetween those of Z2* and Z2R (table [2J figure |6|. 





PARP(82) 


PARP(90) 


a 


Z2 


1.83 


0.28 


0.00 


Z2* 


1.93 


0.23 


0.00 


Z2R 


1.87 


0.23 


4.15 


Z2R' 


1.90 


0.23 


2.67 



Table 2: Result of the 2-step 3-parameter tune to all UE and MB observables. 



3.4 The CMS ridge 

Finally we consider correlation function R{Ar]^A(j))^ comparing results using tunes Z2R and 
Z2R' with those from the original paper (figure ItI (top row)). In the middle row, the results 
for Z2R are shown. It is clear that for high- multiplicity moderate-p^ events (middle, left), the 
long-range near-side ridge is visible, fully in agreement with the CMS results. In the same 
row (middle, centre), one can see that also for moderate- multiplicity events a ridge is visible, 
denoting that perhaps the effect of the modification is too strong. In the bottom row, the 
same plots are shown for tune Z2R'. Here, the effect is not strong enough at high-multiplicity 
(bottom, left), as no near-side ridge is visible, while it is still too strong at moderate-multiplicity 
(bottom, centre), where an unwanted ridge is visible. For high- multiplicity events, including all 
Pt > 0.1 GeV (middle/bottom, right), both tunes show similar effects. There is no near-side 
ridge and some broadening around A77 = is visible, both in agreement with CMS data, but 
there is an unexplained additional peak at (A77, A(/)) = (0,7r). 



4 Conclusions 

We proposed a modification of PythiaG, explaining the ridge effect with multiparton interac- 
tions. The model introduces a correlation between the azimuth of the event planes of individual 
multiparton interactions and the event plane of the hardest interaction. This correlation can 
be naturally explained in a physical picture based on the impact parameter between the pro- 
tons. The two main implications of this modification are the appearance of the near-side ridge 
in high- multiplicity moderate-p^ events and a shift in the activity in the transverse region. 
This latter effect can be counteracted by a re-tune of the pT-cutoff parameters to underlying 
event data. In a slightly broader picture, minimum bias data can be included in the re-tuning. 
Implementing this with the PROFESSOR package, we found tunes Z2R and Z2R'. 
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Nch at yi = 0.9 TeV, \rj\ < 2.0, 60 < \A(p\ < 120 




- Z2* (1.93, 0.23, 0.00) 

- Z2R (1.87, 0.23, 4.15) 

- Z2R' (MB) (1.90,0.23,0.00) 



Leading Track-jet pj- [GeV /c] 
Nch at yj = 0.9 TeV, \t]\ < 2.0, 60 < \A<p\ < 120, pf^" > 3GeV/c 




- CMS UE data 

- Z2 (i.f 

- Z2* (1.93, 0.23, 0.00) 

- Z2R (1.87, 0.23, 4.15) 

- Z2R' (MB) (1.90,0.23,0.00) 

- Z2R' (MB.UE) (1.90, 0.23, 2.67) 





Nch at ^ = 7.0 TeV, \r!\ < 2.0, 60 < \A(j>\ < 120 



- CMS UE data 

- Z2 (1.83, 0.28, 0.00) 

- Z2* (1.93, 0.23, 0.00) 

- Z2R (1.87, 0.23, 4.15) 

- Z2R' (MB) (1.90,0.23,0.00) 



Leading Track-jet p^ [GeV /c] 
Nch aty^= 7.0 TeV, \i^\ < 2.0, 60 < \A(p\ < 120, pf*^" > 3GeV/c 




- CMS UE data 

- Z2 (1.83, 0.28, 0.00) 

- Z2* (1.93, 0.23, 0.00) 

- Z2R (1.87, 0.23, 4.15) 

- Z2R' (MB) (1.90,0.23,0.00) 

- Z2R' (MB.UE) (1.90, 0.23, 2.67) 




Nchat\/i=7.0TeV, \t]\ < 2.4 




- CMS MB data 

- Z2 (1.83, 0.28, 0.00) 

- Z2* (1.93, 0.23, 0.00) 

- Z2R (1.87, 0.23, 4.15) 

- Z2R' (MB) (1.90,0.23,0.00) 

- Z2R' (MB_UE) (1.90, 0.23, 2.67) 




- CMS MB data 

- Z2 (1.83, 0.28, 0.00) 

- Z2* (1.93, 0.23, 0.00) 

- Z2R (1.87, 0.23, 4.15) 

- Z2R' (MB) (1.90,0.23,0.00) 

- Z2R' (MB.UE) (1.90, 0.23, 2.67) 




Figure 6: Performance of the Z2R [red solid] and the Z2R' [green dashed (intermediate), sohd 
(final)] tunes, compared to tunes Z2 [black dashdotted] and Z2* [blue dashed], for select ob- 
servables in the full UE (rows 1-2) and MB (rows 3-4) data set. The four observables used for 
tune Z2R are given in rows 1 and 3. 
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Figure 7: Results for R{Ar]^A(f)): original (top row), Z2R (middle row) and Z2R' (bottom 
row); (left) high multiplicity, moderate pr; (centre) minimum bias, moderate pr; (right) high 
multiplicity, all pr > 0.1 GeV. 
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